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Abstract 

In this technical note we calculate the dynamics of a linear feedback model of pro- 
gression in the cell cycle in the case that the cells are organized into fc = 3 clusters. 
We examine the dynamics in detail for a specific subset of parameters with non-empty 
interior. 

There is an interior fixed point of the Poincarc map defined by the system. This 
fixed point corresponds to a periodic solution with period T in which the three clusters 
exchange positions after time T/3. We call this solution 3-cyclic. In all the parameters 
studied, the fixed point is either: 

• isolated and locally unstable, or, 

• contained in a neutrally stable set of period 3 points. 

In the later case the edges of the neutrally stable set are unstable. This case exists 
if either the three clusters are isolated from each other, or, if they interact in a non- 
essential way. In both cases the orbits of all other interior points are asymptotic to the 
boundary. Thus 3-cyclic solutions are practically unstable in the sense the arbitrarily 
small perturbations may lead to loss of stability and eventual merger of clusters. Since 
the single cluster solution (synchronization) is the only solution that is asymptotically 
stable, it would seem to be the most likely to be observed in application if the feedback 
is similar to the form we propose and is positive. 

January 20, 2013. 



1 Introduction 



This technical note is a supplement to the manuscript [T]. In that work the dynamics of 
cells in a cell cycle were considered when cells are subjected to feedback from other cells 
in the cell cycle and the response to the feedback is also dependent on the cells location 

* Centre de Physique Theorique (UMR 6207 CNRS - Universite Aix-Marseille II - Universite Aix- 
Marseille I - Universite Sud Toulon- Var) CNRS Luminy Case 907, 13288 Marseille CEDEX 9, France 

t Corresponding author. 321 Morton Hall, Mathematics, Ohio University, Athens, OH 45701, USA, 
youngOmath . ohiou . edu 



1 



within the ceU cycle. We refer the reader to that manuscript for the biological motivation 
for this study and for background on the model we study. 

We consider a piecewise affine model for the dynamics of cell populations along the cell 
cycle. Let a population of N cells be organized in k equal clusters {k divides N) labeled 
by a discrete index i G {0, . . . ,k — I}. For the cluster i, the progression along the cycle is 
represented by the periodic variable Xi € M/Z which evolves in time according to the entire 
population status. More specifically, the variables follows the deterministic flow associated 
with a real vector field /, i.e. Xj = fi{, Xi,x, r,s), i = 0, . . . ,k — 1, that is the same for 
each i. Here < s ^ r < 1 are two parameters respectively governing the length of the 
signaling S = [0, s) + Z and of the responsive R = [r, 1) + Z regions. Let a represent the 
fraction of cells in the signaling region S, i.e. 

#{j : Xj e [0, s) mod 1} 
k 

Then we will consider the case that / assumes a simple expression 

dxj _ J 1 if Xi - [xil £ [0,s) 

dt -/^l^'^'^^- I 1 + ^ if ^^_Lx,J G[r,l), 

where [-J denotes the floor function. In short terms, clusters in the responsive region are 
accelerated by a fraction equal to the proportion of cells in the signaling region. In [1] 
we consider systems in which the feedback experienced by cells i? is a general monotone 
increasing or decreasing function of a. 

Relabeling, integer translation of coordinates, and time translation are symmetries of the 
dynamics. Thus, one can assume that all coordinates Xi{0) are initially well-ordered and 
belong to the same unit interval, i.e. we have 

= a;o(0) ^ Xi{0) ^ Xfc_i(0) < 1, i = l,...,k-2. 

The definition of the vector field / implies that the coordinates cannot cross each other 
as time evolves; thus this ordering is preserved under the dynamics. Moreover, the first 
coordinate xq must eventually reach 1 (not later than at time 1), i.e. there exists ^ 1 
such that xo(tR) = 1, and more generally, it must reach any positive integer as time runs. 
Thus the set xq G N defines a Poincare section for the dynamics and the mapping 

(xi(0),X2(0), . . . ,Xfc-l(0)) l-^- (xi(tR),X2(tR), . . . ,Xfc_l(tR)) 

defines the corresponding return map. 

We rely on the following considerations. Starting from t = 0, compute the time ti that 
needs to reach 1 and compute the location of the remaining cells at this date. Define 
F to be this mapping. (We assume first that Xfc_2(ii) < 1 for simplicity, but we will relax 
this soon.) 

F : (xi(0),X2(0), . . . ,2;fc_i(0)) ^ (xo(ti), xi(ti), . . . , Xfc-2(ii)) 
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Figure 1: Illustration of the map F for k = 3. Here F{xi,X2) = {xi,X2) = (xo(ti), 
Notice that xo(ii) = ti by assumption on xo(0). 

Now the time ti +t2 that xi._2 needs to reach 1, together with the population configuration 
at t = ti + t2, follow by applying F to the configuration {xQ(ti),xi{ti), . . . ,Xfc_2(ti)). By 
repeating the argument, the desired return time % is given by = ii + t2 + • • • + ifc and the 
desired return map is F'^. Therefore, to understand the dynamics, one only has to compute 
the first map F. 

1.1 General properties for arbitrary k 

It was noted in |T] that we may regard F as a continuous piecewise affine map of the 
(k — l)-dimensional simplex 



into itself. (Although the boundaries and 1 are identified in the original fiow, in the 
analysis here, we consider them as being distinct points for F.) 

On the edges of the simplex, F has a relatively simple dynamics. Indeed if, initially, all 
coordinates are equal, then they must all reach the boundary 1 simultaneously. In other 
words, on the diagonal (xj = x for all i), we have F{x, . . . ,x) = (ti, !,...,!) where ti 
depends on r,s and x (for x = 0, we have ti = 1 independently of r and s). Moreover, 
starting with x^-i = 1 implies ti = which yields 



whatever the remaining coordinates xi, . . . ,Xk-2 are. As a consequence, the edge 



which is mapped onto {(0,0, x, !,...,!) : x € [0, 1]} and so on, until it reaches the edge 
(0, . . . , 0, x) which is mapped back onto the diagonal (after k iterations). 

A particular orbit on the edges is the A;-periodic orbit passing the vertices, and which 
corresponds to the single cluster of velocity 1 in the original flow, namely 



(0, . . . , 0) ^ (1, . . . , 1) ^ (0, 1, . . . , 1) ^ (0, 0, 1, . . . , 1) ^ . . . ^ (0, . . . , 0, 1) ^ (0, . . . , 0) 



^ Xi ^ X2 ^ . . . ^ X/^._l ^ 1 



F(xi, . . . , Xfc_2, 1) = (0, Xi, . . . , Xfc_2) 



{(x, !,...,!) : XG [0,1]} 



is mapped onto 



{(0,x,l,...,l) : XG [0,1]} 
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From a direct analysis in the original system [T], we know that this orbit must be asymptot- 
ically stable for positive feedback and unstable for negative. Since there are two points of 
a periodic orbit at each boundary of every edge (which are themselves globally /c-periodic 
1-dimensional sets), by the intermediate value theorem, there must be at least one other 
/c-periodic orbit on the edges with coordinates comprised between and 1. Whether this 
orbit is unique might depend on parameters. 

Finally, since the simplex is a convex and compact invariant set under F and the boundary 
cannot contain any fixed point, the Brouwer fixed point theorem implies the existence of a 
fixed point in its interior. From the definition of F this corresponds to a solution satrisfying: 

Xi{ti) = Xi+i(0) for all z = 0, . . . , A; — 2, and = xo(0) mod 1. (1-1) 

We call this type of solution a /c-cyclic solution [1]. 

2 Dynamics for A; = 3 

We will study the dynamics of the linear model only for a limited subset of parameter 
space. It will become clear to the reader that the analysis could be reproduced for the 
entire parameter space, but doing so would require a prohibitive amount of time and space. 

2.1 Computation of the map F 

When /c = 3, F is defined in the triangle ^ xi ^ 2:2 ^ 1 and maps the pair (xi,X2) into 
(ti, xi(ti)), i.e. not only the hitting time ti needs to be computed now, but also the location 
of xi at this instant has to be specified. There are numerous cases to study depending on 
xi and X2. To list them, we primarily use the location of X2. 

• 0^X2<r — s. As for /c = 2, in this case X2 moves with velocity 1, even when 
lying in R. Therefore ti = 1 — X2. Since xi ^ X2, we also have xi ^ r — s and thus 

Xl{ti) = Xl + tl = 1 - X2 + Xl. 

• r — s ^ X2 < r. In this case, X2 enters the responsive region, but only after a while, 
precisely at t = r — X2. Its acceleration depends on the location of xi at this date, 
whether xi belongs to S or not. This alternative is given by the relative sizes of 
r — X2 and s — xi. Considering also the possible cases of acceleration for xi we get the 
following sub-cases. For simplicity, we assume that s < r — s and r -|- |s < 1. 

* ^ xi ^ s and < r — X2 < s — Xi. In this case, xi is still in S when X2 enters 
R. The assumption r -|- |s < 1 implies that both xq and xi get out of S before 
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X2 reaches 1 (existence of the third hne in the expression below). We then have 

X2{t) = 



r + |(t — r + X2) if r — a;2 < t ^ s — xi 



3 

r+|(s — xi — r + X2) + |(t — s + xi) if s — xi < t ^ s 
r + |(s — xi — r + X2) + |xi + t — s if s<t 

which imphes ti = 1 — |x2 + ^ + |(r — s). Furthermore, the assumption s < r — s 
forces xo(i) to be out of 5 when xi{t) enters R (if it does). Then xi(ti) = xi+ti 
in this case. 

* ^ xi ^ s and s — xi^r — X2^sors<xi^r — s and X2 < r. Here, xi(t) is 
out of S when X2{t) enters R. From parameter assumptions we have r + |s < 1 
and thus xq leaves S before X2 reaches 1. Consequently, the evolution of X2 reads 



X2(i) 



r + |(t — r + X2) if r — X2 < t ^ s 



3 

r + |(s — r + X2) + i — s if s^t 



from which we get ti = 1 — 3X2 + and again xi(ti) = xi + ti. 

* r — s < xi ^ X2 and r — s < X2 ^ r. In this case, ti remains as before. However, 
xi now suffers some acceleration and we have 



r + -^[s — r + xi) + t — s if s<t 



Consequently xi(ti) = 1 — |(x2 — xi). 

r ^ X2 and ^ xi ^ s. For X2 larger than r it is useful to consider separately the 
cases xi ^ s and xi > s. In the latter case, we always have xi(ti) = xi + ti but the 
expression of ti depends on X2. We have 3 sub-cases 

* ^ xi ^ s and r^X2 ^ 1 — |,s + ^. Here X2 starts sufficiently near so that 
both xi and xq get out of S before ti 



X2{t) 



X2 + |i if < t < s - xi 

X2 + |(s — xi) + |(t — s + xi) if s — Xl < t ^ s 
X2 + ls-^+t-s if s <t 



which implies ti = 1 — X2 + ^ — |s. 

* ^ xi ^ s and 1 — |s + ^<X2^1 — |(s — xi). In this case, we have ti ^ s 
and hence X2 reaches 1 during the second phase in the previous expression. This 
yields ti = l{l-X2 + f-^). 

* ^ xi ^ s and 1 — |(s — xi) ^ X2 ^ 1. There X2 reaches 1 even before xi 
reaches s. This results in ti = |(1 — X2). 



5 



• s < Xl. When xi > s, the situation for X2 is similar to as for k = 2. This point may 
reach 1 before or after xq leaves S. The first case occurs when r ^ X2 ^ 1 — |s and 
gives ti = 1 — X2 — |. In the second case 1 — |s < .X2 ^ 1, we have ti = |(1 — X2). 

We can now focus on xi{ti). When r ^ X2 < 1 — |s and xi ^ r — s, there cannot be any 
acceleration and thus xi(ti) = xi + ti. Furthermore, when 1 — |s ^ X2 ^ 1, we have ti < s 
and thus xi receives any acceleration only if it is initially larger than r — ti. In practice, 
this leads the following cases. 

• r — s < Xl < r and r ^ X2 < 1 — |s. In this case, s < ti and the evolution of xi is as 
in equation ()2.ip . It results that xi(ti) = 1 — X2 + |xi — |. 

• r — 1(1 — X2) < Xl < r and 1 — |s ^ X2 ^ 1, then only the first phase in (12. ip applies 
and we get again xi(ti) = 1 — X2 + |xi — |. 

• r ^ Xl < X2 and r ^ X2 < 1 — |s. Now, we have 

. , I Xl + it if t ^ s 
xi{t) = < . 

I xi + gS + i — s if s^i 

Thus xi{ti) = Xl + ti + I = 1 — X2 + Xl. 

• r ^ Xl < X2 and 1 — |s < X2 ^ 1. As ti < s, only the first phase applies in the 
previous expression. We get xi(ti) = 1 — X2 + xi. 

All together, when the parameters are such that 2s < r < 1 — |s, the triangle decomposes 
into 13 polygonal subdomains with affine dynamics, see Figure[2j For simplicity, we have la- 
beled these domains by integer numbers following the increasing ordering in the coordinates 
Xl and then in X2. For the sake of clarity, we provide the following table that recapitulates 
the domains and the corresponding expression of F. Since the map F is continuous (but 
not C^), the expressions coincide on the domain boundaries. Accordingly, (when stability 
is not considered) boundaries can be regarded as simultaneously belonging the adjacent 
domains. 
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Label 


Condition on xi 


Condition on xi 


a;o(ii) 


X\(t\) 


1 


.Tl < X2 


^ a;2 < r — s 


1 - X2 


1 — X2+ Xi 


3a 


^ xi ^ s 


r — s^X2<r — s + xi 


1 4X2 1 r — s 
^ 3 3 




2 


^ xi ^ s 


r — s + xi ^X2 <r 


1 5x2-3:1 , 2(r-s) 
3 3 


1 5x2-4x1 1 2(1 — s) 
3^3 


5 


Os^xis^s 


r^a;2<l-f + f 


l-.X2 + f 


l-X2 + ^--f 


10 


^ xi ^ s 


1 ^^71 a;2 < 1 '^-'7^' 


3(1 — X2) 1 xi — s 
4 ' 4 


3(1 — X2) 1 5xi-s 
4 ' 4 


9 


O^xi 


1 '(^7^^ ^ X2 < 1 


3{1-X2) 

5 




3b 


s < xi ^ r — s 


r — s ^ 2:2 < r 


1 4X2 I S 

^ 3 3 




6 


s < Xi ^ r — s 


r a;2 < 1 - ^ 


1 - a;2 - f 


1 - X2 + Xi - 1 


11 


s<xx^r 4 


1 - f 2^2 < 1 


3(1-X2) 

4 


XI + '''-/^> 


4 


r — s < xi ^ 0:2 


r — s ^ a;2 < r 


1 4X2 1 r — s 
^ 3 3 


^ 4(X2-Xi) 


7 


r ~ s < x\ ^ r 


r a;2 1 - ^ 


1 ^ ^2 ~ f 


1 - :E2 + ^ - i 


12 


r ^ 4 ^ <X\^r 


1 - f < X2 ^ 1 


3(1-X2) 

4 


l--2 + ^-f 


8 


r < x\ ^ X2 


r a;2 ^ 1 - ^ 


1 ~ ^2 - f 


1 - a;2 + xi 


13 


r < x\ ^\ 


1 - |s < 12 ^ 1 


3{1-X2) 

4 


1 — a;2 + 2^1 



2.2 The symbolic dynamics 

In order to obtain insights into the dynamics, we study the image of the triangle partition 
by using explicit expressions in the table above. To that goal, we begin with a series of 
elementary observations. 

• We already know from the analysis for arbitrary k that the 3 corner points are mapped 
one into another in a stable period-3 orbit. We also know that there must be another 
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3-periodic orbit on the edges. 

• In addition, there must be a fixed point inside the triangle. 

• F{a) = F(0, r - s) = {1 - r + s,l - r + s) and F{b) = F{r - s,r - s) = {1 - r + s,l). 
So the horizontal segment (a, b) is mapped into a vertical one with abscissa 1 — r + s. 
Depending on the location of r in the interval [2s, 1 — -y] , the quantity 1 — r + s 
varies in the interval [s,l]. Precisely, we have (one or several cases might not apply 
depending on s) 



Conditions on (r, s) 


Location of the segment [-F(a), F{b)] 


2s ^ r ^ 


beyond r 


i±^<r^i + s 


between r — s and r 


^+s<r^l-^ 


between s and r — s 



• Similarly, the segment (d, e, /) is mapped into a vertical one since we have F{d) = 
F{s, r) = (1 - r - §, 1 - r + f ), F{e) = F(r - s, r) = (1 - r - §, 1 - f ) and 
F{f) = F{r,r) = (1 — r — |, 1). The parameter dependent location of 1 — r — | is 
listed in the following table 



Conditions on (r, s) 


Location of the segment [F (d), 


2s ^ r ^ 1 - § 


beyond r 


2 6 ^ ^ 2 ^ 3 


between r — s and r 


2^3^^ 3 


between s and r — s 



In addition, we have F(c) = F(0, r) = (l-r-^,l-r-^). 

• The segment {g, h, k) is mapped into one with abscissa s and in particular F{g) = 
F(0,1 - f ) = {s,s), F{h) = (s,2s), F{i) = (s,r) = d F{j) = (s,r + f ) and 
F{k) = (s, 1) = /. 

• We already know that the segment {l,ra) is mapped onto the vertical axis xi = 0. 
Moreover F{1) = (0, s) and F{m) = (0, r) = c. 

The two previous properties imply that the regions in the upper vertical strip are mapped 
into the following ones (see Figure [2D 

9, 10 1 

11 lU2U3a 

12 5 

13 5U 10U9 

For the remaining regions, the situation depends on parameters and decomposes in various 
cases according to the two previous tables. We consider here the case where all points from 
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a to f fall beyond the vertical line xi = r, i.e. 2s < r ^ ^ — |. From Figure [21 we get the 
following properties 

1,2,3,4 8 U 13 
5,6h^lU3U4U6U7U8 
7h^6U7U8Ul3 
8h^6U7UllUl2Ul3 



2.3 Parameter dependence of fixed points and bifurcations 

The above properties suggest to consider possible fixed points in 6 and in 7. It turns out 
more convenient to study the one in 7 and to follow its bifurcations in parameter space. 

Fixed point in 7: The (unique) fixed point in 7, namely {^-^-j^^, turns out to be 

a source (with associated complex eigenvalues for - the return map under study) which 
exists provided that 

3 - 2s 3 + 8s 



9 9 

At the lower boundary r = ^~g^ of this domain, the first coordinate meets the vertical 
boundary with 8. This suggests that (^-'-fg-^, ^^~^~^) might be continued as a fixed point 
in 8 for r < This is indeed the case. 

• Fixed point in 8: Easy calculations conclude that the coordinates are given by (^-^, ^^)- 
Moreover, the fixed point is neutral (with a double eigenvalue of equal to 1) and 

it exists provided that 

3 - 2s 

2s ^ r ^ 

9 

For r = this point coincides with the fixed point in 7 and we expect a kind of 
pitchfork bifurcation which would create a 3-periodic orbit with code 7 i— )> 7 i— >■ 8 i-^ 7. 

• Period-3 orbit 7 i— >■ 7 i— ?• 8 iH> 7: Calculations show that this (unique) orbit has coor- 
dinates: 

r,l — r ) 1-^ (r, 2r H — ) i-^ (1 — 2r , 1— r jiH^lr, 1— r ) 

3; V , -r 3^ V 3 , 3^ ^ ' 3^ 

and indeed emerges from the fixed point in 7 at r = It persists for all lower 

values of r down to 2s and is expanding with double real eigenvalues. 

• Period-3 orbits 8 1— )• 8 1— > 8 1-^ 8: Topological considerations suggest that other orbits 
should come into play at r = . Indeed, calculations confirm that a two-parameter 
family of 3-periodic orbits lying in 8 also emerges from the fixed point in 7 when r 
crosses downward. All orbits are neutral with double eigenvalue 1 and their 
coordinates have the following expressions 

s s s 

(Xl, X2) ^ (1 - X2 - 3, 1 - X2 + Xi) (X2 - Xi - -, 1 - - - Xi) ^ (Xi, X2) 
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1-r-s/3 



2r+s/3 

r 

Figure 3: Invariant triangle composed of the neutral fixed point and neutral period 3 orbits 
lying in 8. The 3 corners are points of the 3-periodic orbit with code 7 i— t- 7 i— > 8 i— t- 7. 
Notice that the left vertical edge (belongs to 8 and) is part of the boundary between 7 and 
8. 

where xi is arbitrary in the interval (r, X2 — t — |) , X2 is arbitrary in (2r + |, 1 — r — |) 
and r is arbitrary in [2s, ■ The orbits can be viewed as rotating around the fixed 
point (^-^, which is included in the family {i.e. for {xi,X2) = {^-^, ^^)) t^ie 
orbit actually reduces to a fixed point). The closure of the set of orbit coordinates 
forms a triangle whose corners are the coordinates of the orbit 7 i— > 7 i— > 8 i— )■ 7 (see 
Figure [3D. 

At the upper boundary r = of the "7"-fixed point's existence domain, its first coordi- 
nate meets the vertical boundary with 6. We then consider the fixed point in 6. 

• Fixed point in 6: Its coordinates are given by (^^, ^) and this neutral fixed point 
exists iff 

3 + 8s 2(3 -s) 

9 9 

The bifurcation scenario at r = is similar to the one above at the lower boundary 

_ 3-25 

' ~ 9 • 

• Period-3 orbit 7 i-^ 7 i— > 6 i-^ 7: It turns out that this unique orbit emerges from the 
fixed point 7 at r = and exists up to r = (or up to 1 — ^ whichever 
is smaller). It is expanding with double real eigenvalues. At r = both the 
component in 6 and its successor cross the horizontal line X2 = r and the second 
component in 7 cross the upper domain boundary X2 = 1 — ^ . 

However, the value r = does not involve any existence condition related to the 
fixed point in 6. As shown below, the reason is that the periodic orbit can be continued 
up to ^^^"'^^ provided that we consider the appropriate sequence of symbols. 
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Period-3 orbit 6 i-^ 3 i— ?• 3 iH> 6: This orbit continues the orbit 7 i— ?> 7 i-^ 6 i— )> 7 when 
r > Indeed, it exists provided that 

3 + 2s 2(3 -s) 

6 9 

It is expanding with 2 real eigenvalues (the same as those associated with 7 i— >■ 7 i-^ 
6 I—)- 7) and coincide with 7i— 7>7i— >6i— )-7atr = '^^^^ . 

Period-3 orbits 6 i— )■ 6 i— ?• 6 iH> 6: Similarly to as before, a triangle of neutral 3-periodic 
orbits is created from the fixed point 7 at r = • The coordinates have the following 
expression 

{xi,X2) 1-^ (1 - X2 - |, 1 - X2 + - |) 1-^ (2:2 - Xi, 1 - j;i - |) (xi, X2) 

and these orbits exist in the same parameter domain as the fixed point in 6. When 
^ £ [^^^5 ^^^" j the restrictions on coordinates are 

-, o 5s , ^ 2s 

1 — 2r -\ < xi < r — s and 1 — r -\ < X2 < xi + r — s 

3 3 



and they are 



when r € 



s s s 

2r -\ l<xi<l — r and r < 2:2 < xi + 1 — r 

3 3 3 

3+2s 2{3-s) 



6 ' 9 



At the boundary r = the expanding periodic orbit and the family of neutral orbits 

meet with the neutral fixed point "6" in an inverse pitchfork bifurcation to create the fixed 
point in 3. 

• Fixed point in 3: The fixed point {^^jf^, Hl^zlzfl) g 3 is expanding with double real 
eigenvalues and exists in the domain 

2(3 -s) 2 
9 3 

At the upper boundary r = | + s of its existence domain, the fixed point "3" meets the 
region 1. 

• Fixed point in 1: The fixed point (5, |) is neutral with double real eigenvalue 1 and 
exists in the domain 

2 

- + s < r < 1 
3 
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0.0 0.1 0.2 0.3 0.4 0.5 

S 



Figure 4: A; = 3. Parameters domains of existence of fixed points and related periodic 
orbits analyzed in this section. In regions (2) and (4) the unique fixed point is unstable. In 
regions (1), (3) and (5) the fixed point is neutral and sits inside a triangle of neutral period 
3 points. 

• Period-3 orbit 3 i— >■ 3 i— ?• 1 iH> 3: Same parameter domain. Expression of coordinates 
{2{r - s) -l,r - s) {l-r + s,r-s) (1 -r + s, 2(1 -r + s)) i-^ (2(r - s) - 1, r - s) 
Expanding with double real eigenvalue. Emerges from the fixed point 3. 

• Period-3 orbit 1 i-^ 1 i— t- 1 i-^ 1: Same parameter domain. Triangle of neutral 3-periodic 
orbits 

{XI,X2) H> (1 - X2, 1 - X2 + 1-^ (X2 - Xl, 1 - Xi) (xi, X2) 

where xi is arbitrary in (1 — r + s, 2(r — s) — 1) and X2 S (1 — r + s + xi,r — s). As 
before, the corners coincides with the period-3 orbit 3 i— > 3 i— )■ 1 i-^ 3. 

To recapitulate, the return map possesses the following orbits depending on parameters 
in the considered regions (see Figure H]) 

(1) 2s ^ r ^ three sources (two of them belong to 7, the other one lies in 8) and a 
two-parameter family of neutral fixed points lying inside the triangle whose corners 
are the 3 sources. 

(2) ^ < r ^ One source in 7. 

(3) ^ < r ^ Similarly to as in (1); namely three sources (two in 7, one 6 if 
r ^ ^"^^^ , and two sources belong to 3 and one is in 6 for r ^ ^-^) and a triangle of 
neutral fixed points in 6. 
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(4) < r ^ I + s. Similarly to as in (2); namely a single source in 1. 

(5) I + s < r < 1 — Similarly to as in (1); namely three sources (two lie in 3 and 1 
lies in 1) and a triangle of neutral fixed points in 1. 

Note that in case (5) we have r — s > | and this corresponds to the case that the three 
clusters in the cyclic solution (which has initial conditions (0, 1/3, 2/3) do not interact, i.e. 
xq leaves S before X2 enters R and no feedback is experienced. 

In the cases (1) and (3) the clusters in the cyclic solution experience feedback, but in a 
non-essential way, by which we mean that small perturbations do not lead either to further 
separation or contraction between the clusters. For instance in case (3) X2 begins in R and 
xi begins in between S and R and xq leaves S before either of these states changes. In case 
(1) both xi and X2 begin in R, but xq leaves 5* before either leaves R. 

A python script that will produce movies for the dynamics for arbitrary < r < s < 1 can 
be found at: http://oak.cats.ohiou.edu/~rb301008/research.html. 

2.4 Additional orbits 

Besides orbits bifurcating with fixed points, based on numerics and on properties of the 
symbolic dynamics, additional orbits of F exist depending on parameters. Their existence 
domains do not coincide with those listed above. 

• The most important orbits are the one-parameter family of period-3 neutral-stable 
orbits lying on the edges (and with code li-^-li-^lli-^l) 

(0, X2) ^ {I - X2,l - X2) {X2, 1) 

These orbits exist for arbitrary X2 G {1 — r + s,r — s) provided that ^ r < 1 — ^• 
They are stable with respect to transverse perturbations. The coordinates form 3 
open intervals, each included in one the edges. Moreover, the interval boundaries 
respectively form a period-3 hyperbolic orbit (code 2i-^li-^lli— )-2) and a period-3 
neutral- unstable orbit (code li— )-4i-^lli— >1) which exist in the same parameter 
domains and which merge for r = ^ + s. Orbits in the family correspond to 2 clusters 
solutions in the original system with one cluster composed of two clusters (and the 
two clusters being certainly isolated one from each other). 

• Interestingly, the bifurcation at r = ^ -|- s generates two distinct periodic orbits with 
identical code (2 1— > 4 1— )• 11 1— > 2). Both orbits are neutral-unstable and the first has 
coordinates 

,3-6(r-s) 3-2(r-s) 3 - 2(r - s) 

( \ -, \ -) ^{r-s, ^{r-s, 2{r - s)) 
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with the first point being in 2. It exists under the condition ^ + |^r^^ + s. 

The other solution is actually a one-parameter family of periodic orbits (which forms 
a segment and) whose component in 11 is (^^^^—j^-^^,X2^ where X2 can be chosen 
arbitrary with the condition 

3 ( 3 — 6(r — s) 6r + 5s — 3 

— (1 — X2) ^ mm ' 



The family exists under the condition i — ^^r^^ + s. 

The bifurcations taking place at r = ^ + | and f = ~ ^ ^-I'e unclear. The fact that the 
one-parameter family reaches the boundary with 5 in the former case suggests to investigate 
orbits with code 5i— )-4i— )-lli— 7-5. Surprisingly, the analysis concludes that such an orbit 
exits only if r = ^ + §, a value that is unrelated to the one-parameter family existence 
condition. 

Two additional orbits have been found on the edge when 2s < r ^ ^ -|- One is hyperbolic 
with code 2i-^'8i-)-lli-)-2 and coordinates 

(0, r - ^) ^ (1 - r + ^, 1 - r + ^) ^ (r - |, 1) ^ (0, r - ^). 

The other one is neutral-unstable with code 2i-t'4i— )-lli— )-2 and coordinates 

(0, ^-l^tzll) ^ (1 _ ^+^-\ 1 _ ^-±1-^) ^ (^±^^, 1) ^ (0, ^ + ^1^-^^ 

The two orbits merge for r = ^ -|- ^ in a seemingly saddle-node bifurcation. 



3 Discussion 



In this note we have investigated only a portion of the possible parameter space < s < r < 
1. It should be clear to the reader by this point that further investigations into other subsets 
of the parameters in this fashion are possible, but perhaps prohibitively time-consuming. 
We see that such studies are likely to also prove unprofitable, since numerical simulations 
show that no other types of dynamics occur other than those described here. One can 
download a python script that to investigate the dynamics for arbitrary 0<r<s<lat: 
http : //oak. cats . ohiou. edu/~rb301008/research.html. 

The dynamics we have observed for these parameter sets closely resembles the dynamics of 
k = 2 cluster systems analyzed in [1]. The k = 3 fixed point of F with positive feedback, 
like that for A; = 2 is either: 



• isolated and locally unstable, or, 

• contained in a neutrally stable set of period k points. 
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In the later case the edges of the neutrally stable set are unstable. This case exists if either 
the three clusters are isolated from each other, or, if they interact in a non-essential way. 
In both cases the orbits of all other interior points are asymptotic to the boundary. Thus k 
cyclic solutions for either k = 2 or k = 3 are practically unstable in the sense the arbitrarily 
small perturbations may lead to loss of stability and eventual merger of clusters. Since the 
single cluster solution (synchronization) is the only solution that is asymptotically stable, 
it would seem to be the most likely to be observed in application if the feedback is similar 
to the form we propose and is positive. 
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